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■ Two twisted doublets, one containing the up and down quarks and the other con- 

i taining the strange quark with an SU(2)-flavor partner, are used for studies in the 

P3 ' meson sector. The relevant chiral perturbation theory is presented, and quenched 
QCD simulations (where the partner of the strange quark is not active) are per- 

o . 

CN . formed. Pseudoscalar meson masses and decay constants are computed; the vector 

CN \ and scalar mesons are also discussed. A comparison is made to the case of an un- 



m 
o 

o 

o 



X 



twisted strange quark, and some effects due to quenching, discretization, and the 
definition of maximal twist are explored. 

I. INTRODUCTION 



Oh. Twisted mass lattice QCD (tmLQCD) is a variation on the Wilson action — essentially 



a chiral rotation of quark flavor doublets, acting on quark mass terms relative to Wilson 
terms in the action — whichproduces two desirable features: the removal of unphysical zero 



modes in quark propagator_s 
spacing) at maximal twist 



:s i| 

m 



l| and the elimination of 0(a) artifacts (where a denotes lattice 
A number of numerical simulations have been performed for 
both quenched and dynamical tmLQCD (for a recent review, see Ref. Q]). As well, the chiral 
perturbation theory for tmLQCD (tm^PT) has been developed. It differs from continuum 
xPT by discretization effects and is required for the extrapolation of tmLQCD data. The 
effective theory has also played a vital role in understanding various aspects of tmLQCD such 



as 0{a) improvement, t 



le phase diagram, and the relationships between various definitions 



of maximal twist I 



IG 



With an interest in the phenomenology of hadrons built of m, d and s quarks, our goal 
in this paper is to explore the usefulness of tmLQCD and tm^PT as applied to strange 
hadrons. There is no unique way to introduce the s quark into the calculation. The method 
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used here, to consider a pair of quark doublets {u, d) and ("c", s), is similar to the proposal 



of Pena et al. 



12| |. For the quenched simulations considered here the partner of the s quark 



does not play an active role and should not be thought of as the physical charm quark. 
In this work no mass splitting is introduced within either doublet. The construction of 
the corresponding tmyPT formalism is a straightforward generalization of the published 
one-doublet formalism [7|, |9|. 

As noted above, applying a relative chiral twist has some valuable consequences but there 
are also some less desirable features that have to be dealt with. The tmLQCD action violates 
parity conservation so, in general, correlation functions contain contributions from states of 
both parities. Parity mixing can complicate, in particular, the extraction of matrix elements 
but this can be ameliorated by appropriate tuning of the twist angles. The tmLQCD action 
also breaks the flavor symmetry. For the version of tmLQCD used in this work the members 
of the quark doublets are degenerate in mass but are distinguished by having opposite chiral 
twists. This can lead to mass splittings within hadron isospin multiplets. It will be seen that 
charged and neutral kaons can acquire a mass-squared splitting which is roughly proportional 
to a?. 

To optimize the elimination of Oia] lattice discretization errors one has to tune the 

Q 

chiral twist angles There is not a unique way to achieve maximal twist as has been 
discussed from the point of view of both effective theoryfl Q and simulationQ, 

A standard method for defining maximal twist uses a tuning procedure which 
involves the correlators of the first two isospin components of vector and axial operators 
with the pseudoscalar density 0, [3, ll^- Using two variations of this method, we examine 
the mixing between the third isospin components of scalar and pseudoscalar correlators. 
Ideally one would like to have a tuning to maximal twist which would banish the physical 
pseudoscalar meson from appearing in the wrong parity correlator; the scalar meson with 
its quenched r^'vr^ contribution would similarly be banished from the other parity correlator. 
This is seen not to happen in our simulations. The mixings observed in actual simulations 
represent higher order discretization effects which differentiate between vector-axial tuning 
and scalar-pseudoscalar tuning. 

In this work, we mainly use maximal twist in the doublet containing the strange quark as 
well as in the (u, d) doublet. An alternative procedure would be to set the twist angle for the 
strange quark to zero or equivalently for the quenched theory to consider a twisted (u, d) 
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doublet and a flavor- singlet Wilson strange quark. The latter approach may be a viable 
one for doing full dynamical simulations. The twisted and untwisted strange quark actions 
lead to different patterns of parity mixing and flavour symmetry breaking at non-vanishing 
lattice spacing. We present some results obtained with an untwisted strange quark action 
to illustrate some of these differences. 

At this point it is worth noting that there exist even other approaches for dealing with 
the strange quark. The proposal of Frezzotti and Rossi in Ref. [17[ allows for a nondegen- 
erate doublet in a way which is suitable for dynamical simulations^^. In the limit where 
quark masses are degenerate within each doublet, it is equivalent to the scheme used in this 
paper. However, for nondegenerate quarks, twist and quark-mass splitting are associated 
with different flavor transformation generators. The fermion action contains terms which 
mix flavors so that flavor symmetry breaking effects would be more complicated to deal with 
in simulations and in the effective theory than for the tmLQCD action considered in this 
work. A further example is Ref. where options for tmLQCD chosen to facilitate the 
calculation of the so-called kaon bag parameter are discussed. 

In addition to meson masses, the pseudoscalar meson decay constants are also consid- 
ered. With quark masses flxed by physical meson masses, the decay constants fj^ and fx 
become absolute predictions, and are shown to compare favorably with previous quenched 
simulations using other actions. All results are consistent with tm^PT. 

The remainder of the article is organized as follows. Section ITTl deflnes the effective chiral 
Lagrangian with two twisted flavor doublets, and Sec. Illll uses that Lagrangian to derive ex- 
pressions involving the pseudoscalar masses and decay constants. Section lTVl presents the tm- 
LQCD action and explains the parameter choices for our numerical simulations, then Sec. |3 
discusses results obtained for the pseudoscalar and vector mesons. Scalar-pseudoscalar mix- 
ing is studied in Sec. IVH and a direct comparison to kaons built from untwisted strange 
quarks is given in Sec. I VI II Section IVIIII contains the conclusions of our work. Details of 
currents and densities in tm^PT are collected in the Appendix. 
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II. THE EFFECTIVE CHIRAL LAGRANGIAN 

To build four-flavor tmxPT, we begin from tmLQCD with two quark doublets, 

I I , i^h= \\ , (1) 



referred to as the light and heavy doublets respectively. Note that the choice of flavor labels 
is a convention; in Ref. [l^ for example, a different choice is made. In this work each doublet 
is taken to be degenerate, so the c quark, which is not active in any of our quenched tmLQCD 
simulations, should not be viewed as the physical charm quark. Pena et a/.|l^ discuss the 
extension of this case to the case of a nondegenerate doublet where the quark-mass splitting 
is aligned with the twist, preserving the favorable feature of no flavor mixing. The fermion 
determinant does not remain real under this generalization so this would not lead a suitable 
action for nonquenched simulations. However, this action may still be useful for valence 
quarks in a mixed action scenario as discussed, for example, in the context of tmLQCD in 
Ref. j2Q|. 

In the so-called "twisted basis" , the two-doublet lattice action is simply a block- 
diagonal version of two copies of the one-doublet theory (the form of which can be found in 
Refs. HQ): 

S'f = <^'Y. ^(^) E ^^^K + - ^ E V.Vm + ^0 + nsMo] , (2) 



where and V* are the usual covariant forward and backward lattice derivatives respec- 
tively, and 

* = , mo = ' , Mo = , (3) 




mh,ol2l \ Hh,or3 

with In the n-by-n identity matrix. The matrix acts in (two-)flavor space and is normal- 
ized so that r| = I2. The parameters nip Q and fip^o are the normal bare and twisted masses 
respectively, with p = I, h. 

Applying the now familiar two-step procedure of Ref. ? an effective chiral Lagrangian 
describing the low energy physics of tmLQCD with two degenerate quark doublets can be 
built as a straightforward generalization of the one-doublet case detailed in Refs. 2,y,|2],Q|- 



From a similar analysis described in Ref . [3| , the form of the effective continuum Lagrangian 
at the quark level is found to be identical to that in the one-doublet case: 

£eff = + ^{P+ m + ii)^! + hsw d ma^^.F^,^ + 0{a^) , (4) 

where Cg is the continuum gluon Lagrangian, and the physical normal and twisted mass 
parameters, m and /x, are defined analogously as in the one-doublet case: 

I • - |_|'^"..'(m,,o-ft.,) " \, (6) 

Zm,himhn - rrich] 




Z^^hfJ'h,o'r3 J \ Zp^f^fih,oT3 

with Zp ; and Zp^h being the matching factors for the pseudoscalar density. The quantities 
fhcj and fhc^h are the critical masses, aside from an 0{a) shift (see Refs. |7|, |9|, |2l|, |22 1 and 
discussions below). Lattice symmetries forbid additive renormalization of fii^o and fih,o- As 
an aside, we note that symmetries also cause the ultraviolet divergent parts of i and h 
to be identical. One can choose a definition of maximal twist (it will be called method (ii) 
in Sec. IIV|1 for which fhcj = fhc^h, but here we do not restrict the discussion to that special 
case. 

Working to NLO in the power counting scheme, 

~ m,, ~ /i, ~ /i/j ~ ~ oAqc-p, , (7) 
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the effective chiral Lagrangian found from matching reads 

£^ = ^Tr(D^SD^St) - ^Tr(xtS + St^) - ^Tr(its + Sti) 

+ L4Tr(D^StD^S)Tr(xtS + St^) + LsTr [(Z}^SD^St)(xSt + Sx^)] 

- Le [Tr(xtS + St^)] ' - Lj [Tr (x^S - St^)] ' - LgTr [(x^S + Stx)^] 
+ iL9Tr( VD^SD.St + + LioTr( 

+ Vr4Tr(D^StD^S)Tr(its + sU) + H^gTr + si^)] 

- WeTiix^i: + Stx)Tr(its + sU) - W;^[Tt{A^^ + SU)]^ 

- WrTiix^i: - S"fx)Tr(i"fS - SU) - Wj[Tt{A^^ - SU)]^ 

+ w,oTt{d^A^d^i: + D^stD^i) 

+ H{Ii{L^,L^, + - i/2Tr(x^x) - H';ii{A\ + x^i) - i/3Tr(iti) , (8) 

where the S field is now SU{A) matrix- valued, and transforms under the chiral group 
5'f/(4)r X SU{A)u. Note that has basically the same form as that in the SU{2) 



22| . except that terms linearly dependent under SU{2) are no longer so 



theory [ 
under SU{A). 

The quantities x ^ind A are spurious for the quark masses and discretization errors, 
respectivelvj23|. At the end of the analysis they are set to the constant values 

X — >2Bo{m + ifi), A — ^21^0^14, (9) 

where Bo and Wq are unspecified constants having dimensions [mass] and [mass^] respec- 
tively. Notice that A involves a single flavour- independent Pauli term for both doublets. 

The discretization effect due to the Pauli term, i.e. the term containing bsw iii Eq. 
can be included non-perturbatively as in Ref. by using the shifted spurion x' = X + ^) 
which corresponds at the quark level to a redefinition of the normal quark mass from m to 

mp = mp + aWo/Bo p = l, h. (10) 

This shift in turn corresponds to an 0{a) correction to the critical mass, so that it becomes 

rricp = Zm,pfhc^p - aWo/B^ p = l, h. (11) 
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In terms of x', the chiral Lagrangian, Eq. Q, can be written as: 

= ^Tr(D^SZ}^St) - ^Tr(x'tS + St^') 

- Li [Tr(Z}^SD^St)] ' - LsTr (D,,SD.St)Tr(D,,SD.St) - LgTr [(D^SZ^^St)^] 
+ L4Tr(D^StZ}^S)Tr(x'tS + St^') + L,Tt [(D^SZ}^St)(x'St + Sx'^)] 

- L6[Tr(x'tS + St^')]' - i^7[Tr(x'tS - St^')]' - i^8Tr[(x'tS + J^^'f] 
+ iL9Tr( VD^SD.St + R^,D^J:^D,^) 

+ W^4Tr(D^StD^S)Tr(its + sU) + H^gTr + si^)] 

- ^y6Tr(x'^S + Stx')Tr(its + sU) - W?^ [Tr(its + sU)]^ 

- WrTiix'^^ - S"fx')Tr(i"^S - SU) - W^[Tr{A^J: - SU)]^ 

- W^8Tr[(x'^S + + SU)] - W^Tt[{A^^ + Sti)^] 

- H^Trix'^x') - H'.TriA^x' + x'^i) , (12) 

where terms that lead only to contact terms in correlation functions and are hence not 
needed below, viz. the Lio, Hi and terms, have been dropped. We have also introduced 
useful combinations 

Wi = Wi-Li, ^ = 4, 5 

Wj = Wj - 2Lj , Wj = W- - Wj + 2Lj , J = 6, 7, 8 

H', = H',-H,. (13) 

As noted in Ref. ^1, the Wiq term is redundant. It can be transformed away by the 
change of variables 

5S = ^^(Sits-i). (14) 



This transforms the Wiq term into the W^, Ws, and H2 terms with their coefficients shifted 
to W5 + Wio, Ws + Wio/2, and H2 — Wiq, respectively. All physical quantities must depend 
then only on these combinations and not on Wio, W5, Ws, and H'2 separately. We have kept 
the Wio term because it provides a useful diagnostic in tm^PT calculations. 
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III. CHIRAL PERTURBATION THEORY FOR GENERIC SMALL MASSES 

In this section, we work out the consequences of the SU(4) effective chiral Lagrangian, 
Eq. ()12|) . which generahzes the results of the SU{2) theory of Ref. Q], and our focus will be 
on the masses and decay constants of kaons: pseudoscalar mesons that involve both flavor 
doublets. We work in the "generic small mass" regime defined by 

Al:^M',>Ml>2Woa, (15) 

where = Anf and 



M; = 25o^m;2 + /i2, p = l,h. (16) 

Note that has dimension [mass^] . In the analysis of our numerical data, quenching effects 
will also be considered (see Eqs. (jlHjl - P^ ). 

A. The vacuum 

At LO the discretized Lagrangian retains its continuum form, so the LO vacuum expec- 
tation value (VEV) of S is that which cancels out the twists in the shifted mass matrix: 

(0|E|0>„.E..^''>^''-'»^'> ° I, (17) 

where Up^o are defined by 

Cp,o = cos(u;p,o) = 2Bom'/Mp , Sp,o = sin(u;p,o) = 25o/ip/Mp , p = l, h. (18) 

This provides one definition for the twist angles. 

At NLO, the VEV of S is realigned by a small amount from Eq. Defining 

/nivln\ -V _ ( exp(z^,,„r3) \ ^ .-,o^ 

{OmO}NLO = ^m = \ , UJp,m = ^p,0 + , P = I, h , (19) 

y exp{iuJh,mT3) J 

the shifts from LO, ei^h, are found from minimizing the potential to be 

ep = - ^^^p'' S^2W, (M; + MO/M; + Ws + ^ [K (Q,o + c.,o) + K Cp,o] } . (20) 

Expanding about the VEV as in Ref. j^, the physical pion fields are defined by 

exp{iLUi m,T3/2) \ 

^ = ^m^ph^m, = I ' , T^ph = exp{i^ / f) , (21) 

exp{iuJh,mr3/2) I 
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where $ has the representation 



1 1 

i=l 



V2 V2' 



12 

V K- k' D- -;|r^8 + ^W 

(22) 

Our choice of the fifteen generators of S'f/(4), Ai, . . . , A15, differ shghtly from the conven- 
tional ones. Here, all off-diagonal generators and the diagonal A3 are the same as the 
conventional ones, but we choose the rest of the diagonal generators to be such that the 
(3, 3) and (4, 4) entries of the diagonal Ag and A15 are interchanged with respect to the 
conventional ones^. This maintains consistency with the ordering of quark fields (m, ci, c, 
s), used throughout this work^, and allows for the standard meson naming convention. In 
particular, we have for the fiavor-diagonal components of $: 

|vr°) = l^[\uu) - \dd)^ , (23) 
Ir^s) = [\uu) + \dd) - 2\ss)^ , (24) 
Ir/is) = + + |ss) - 3|cc) ) . (25) 



12 

Note that the c quark is mass-degenerate with the s quark, so pairs of i^'s and i^'s related 
by the interchange of c and s quarks are mass-degenerate at LO in the chiral expansion. By 
inserting the above expansion of S into the chiral Lagrangian, Eq. ()12j) . the Feynman rules 
for the SU{A) theory can be straightforwardly obtained. 

B. Defining the twist angle 

In the continuum, the twist angle for each degenerate doublet can be defined unambigu- 
ously by 

cUp = tan~^(/ip/mp) , p = l,h. (26) 



^ Note that once the choice is made for Ag, A15 is fixed by the normalization condition Tr(AiAj) — 2Sij. 
^ Recall that u and c have a positive twist whereas d and s have a negative twist. 
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Given the twist angles, the off-diagonal S'f/(4) components of physical currents and densities 
are related to their counterparts in the twisted basis by (with a "hat" denoting the physical 
basis) 

\> = cos(A«a;)\/; + r]at sin(A''^)4 , 
A; = cos(A'^a;)A;: + T],, sin(A'^a;)V;^ , 

= cos{i:''uj)S^ + i sin(S'^cj)P'^ , 
f"^ = cos(S"cu)P'^ + i sin(S"cj)^" , 

a,beR\D, {1,...,15}, 2) = {3, 8, 15}, (27) 



where 



and 



, ±1, for b = a±l 
Vab = { (28) 
0, otherwise 



A'^UJ = -{uJi^ - UjJ , = -(^i, +^j-J , ta<Ja- (29) 



1 1 

The indices ia and ja are the row numbers of the non-zero entries of the SU (4) generator 
Aa that defines the particular flavor current or density, and 



UJi = -UJ2 = ^^l, ^3 = -1^4 = ■ (30) 

For a thorough discussion of currents and densities, see App. |X| We note now that Eq. 
has the form of the inverse transformation of the LO operator, Eq. (jAljl . except that here 
the twist angles are ujp not ujp^m- 

On a lattice, discretization errors mean that different definitions for the twist angles will 
lead to observables that differ by 0{a). In this work, we will define ujp non-perturbatively 
as m Refs. Q,Q,Q,Q by enforcing the absence of parity breaking in the physical basis. 
In particular, we will enforce 

{V^{x)P''{y)) = 0, a,bE^\^. (31) 

From the definitions in Eq. ()27|1 . this condition gives 

{v^{x)P\y)) ^ {v'A^)P'\y)) ,,,, 

tancu; = , , , , , tancu/i = , , ^-,0/ n> • (32) 

(Al(x)Pl(l/)) (A13(x)P13(^)) V ) 

The results for ujp depend on the distance \x — y\ at 0(a). We will enforce the condition in 
Eq. ()31|) at long distance where the single-meson contribution dominates. 
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Evaluating Eq. (j22I) using the results in App. |X1 we find at LO ujp — ujp^yyi — cjp Q, Since 
only and = P", a G {1, 13}, couple to the single-meson state. At NLO, the sole 

non-trivial contributions surviving in the ratios of Eq. (jH^ come from the W\q term, just as 
in the one-doublet theory, and as in Ref. we find 

sinwp,^ ^W^a 
tanwp = ■, = —p^Wio , p = l, h. (33) 

C. Kaon masses and decay constants 

With the Feynman rules in hand, and the twist angles defined, we have now all that is 
needed to calculate pseudoscalar meson masses and matrix elements. At NLO, we find that 
the mass of the neutral kaon is given by 



m 



= M' + ^{4XK]-JKJ} 
+ 4|m'2 (SLe - 4L4 + 2L.-L 



where 



_ _ _ _ 

+M'a (q + Ch){mQ - + - w^) 

+21^2 a2 [(4W?^ + w'^){c, + cnf - K {si + Sh)'] } , (34) 



M' = (M; + M;)/2, (35) 
ml = (M; + 2M0/3, (36) 



m: 



,2 

'Vl5 



(M/ + 5M0/6, (37) 



are respectively the LO K, r/g and 7715 squared masses, and Cp = cosup, Sp = sintUp, which 
we can use instead of coscjp^ and sin rn . respectively at the order we work. Note that the 
usual continuum one-loop contribution 25] appears in Eq. (j34p . 

with Afi being the renormalization scale. Note also that, as required, the neutral kaon mass 
depends only on the combination 2Ws — rather than on Ws and separately, and that 
the kaon mass is automatically 0{a) improved at maximal twist (cp = 0). 
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Flavor breaking in the kaon masses at NLO is given solely by the analytic contribution, 
just as in the SU{2) theory jol, and it reads 

9 9 64^^^9 9 ~, ~, 64:Wn LLl lih 

<o - = --W^ a'siSnW^ = -W^ - ^ , (39) 

where the second equality is derived using the fact that we can replace Sp by Sp^ at the 
order we are working. 

Before calculating Eq. explicitly, one might have anticipated an expression that was 
quadruply suppressed in our power counting scheme, Eq. 0, due to the requirement of 
being O(a^) and 0{fiifih)- However, Eq. shows that the /i^/i/j dependence enters as a 
ratio with M[M'i^ and the squared kaon mass difference is therefore nonzero already at NLO. 

With the physical axial current defined in App. El the decay constant to NLO is 
determined to be 

Ik = /|l + -^[(4^4 + L,)M' + W,a{ci + Ch){m^ + W, + W^,)] 

- \x[ml] - h[ml] - 2JK.] - ix[m^J } , (40) 

where = M/, m|- = mj-, = M' and m|)__ = are the LO expressions for the pion, 
kaon and Dg squared masses respectively. Our tm^PT conventions are such that /^r ~ 93 
MeV. Note that the one-loop contributions from the pions and the r^g are the same as in the 
continuum SU{3) theory 25|. Flavor breaking effects enter first at O(a^), which is NNLO 
for decay constants. The above result shows that the decay constant depends only on the 
combination W5 + Wiq, and is automatically 0(a) improved at maximal twist. 

IV. SIMULATION DETAILS 

We have performed quenched simulations using the action of Eq. (0)- The ensembles 
computed in Ref. [l^, containing 300 gauge configurations each at /3 = 5.85 and jS = 6.0, 
W .Went, "een extended to .L.ude .00 co„fl.u.at,o,. eachfl. 

An additional 

ensemble at /3 = 6.2 has also been generated, again using a pseudo-heatbath algorithm 
which acts on all SU (2) subgroups. Quark propagators are obtained from a 1-norm quasi- 
minimal residual algorithm, with periodic boundary conditions in all directions. Simulation 
parameters are collected in Table HI 
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To remove 0(a) errors, simulations must be done at maximal twist, and to that end we 



employ Eq. (|32|) implemented, following Refs. 



id by 



tancj/ = -— — 7^— j — ^5 , (41) 

where O""^^^ = (O" + iO^)/\/2. This method of tuning separately at each twisted quark 
mass is only one possible choice. A nice summary of the situation in t he g eneric small mass 
regime (Eq. (jl5p ) that we are using has been given by Sharpe in Ref. [ll[. In the language 
of that article, we are using method (i). A variant, called method (ii), extrapolates the 
results of method (i) to the point of vanishing twisted mass, and then uses that definition of 
maximal twist for all values of yUp o (p = l,h). Another option, called method (iii), relies on 
twistings in the scalar-pseudoscalar sector rather than the vector-axial sector. Method (iv) 
assumes that maximal twist can be sufficiently well defined by simply holding the hopping 
parameter fixed at its critical value from the untwisted Wilson theory. 

As sketched in Fig. 1 of Ref. methods (i), (ii) and (iii) are all acceptable non- 
perturbative definitions of maximal twist, and are superior to method (iv). Though most 
of our simulations use method (i), we will make frequent comparisons with results from 

nn 

the xLF collaboration |22l, |2S| using method (ii). We will also present our own results from 
method (ii) when discussing aspects of scalar correlators in Sec. IVII 

Throughout this work, only local operators are used, and error bars reported in graphs 
and tables are statistical only. All statistical uncertainties are obtained from the bootstrap 
method with replacement, where the number of bootstrap ensembles is three times the 
number of data points in the original ensemble. Masses and decay constants are obtained 
from unconstrained three-state fits to correlators, using all time steps except the source. In 
the following sections, each discussion includes references to the relevant figures, but we note 
here that the corresponding numerical values are collected in Table |n] 



V. PSEUDOSCALAR AND VECTOR MESONS 



Masses of the charged and neutral kaons, i.e. the ground state pseudoscalar mesons 
containing one s (anti)quark from the heavy doublet and one u or d (anti)quark from the 
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light doublet, are plotted in Fig. [H Both twisted masses, a/x; o and a^^ Q, take on all values 
from Table |l] and the normal mass term is tuned accordingly (also shown in Table HJ to 
achieve maximal twist at each particular twisted mass value. 

These data a.e expected to be con.ste.t w.tK a ,ue,.hed ^..^MUU of E,, », 
which at maximal twist reads 



M' 



5, 



quench 



In 



Ml 



+ 



ML 



SM'"^ 

+ ^:t^(8L6 - 4L4 + 2L8 



(47r/)2 M;-M/ M'J\ 
64 



U) - 



f 



M' 



64 
64 



w^wi, + o(m;2) + o(M' Vnch) 



M' + 0(M;2) + 0(MVnch) 



P 



L h , 



(43) 
(44) 

(45) 



where ^quench is the standard coefficient for the quenched logarithm. As Fig. [H indicates, lin- 
ear fits in M' to the data at each /3 value yield excellent results, so 0{Mp) and 0(M' 5quench) 
effects are largely unnecessary. Notice in particular that our linear fit for the charged kaon 
did lead to a nonzero (but small) residual kaon mass at M' = 0, thus reminding us that 
corrections to the linear form are important for such details. We have verified, for example, 
that the function AM' + BM' InM' yields an equally excellent fit to our data, and of course 
it enforces the absence of any residual mass at M' = 0. 

Fig. n also reveals differences between our results with method (i) and results from the 

n 

xLF collaboration [22[ using equal quark and antiquark masses with method (ii). For charged 
mesons, the method (i) mass difference is smaller than the method (ii) difference, and as 
expected from tm^PT the chiral limits appear to be very similar. For neutral mesons, the 
chiral limits appear to be somewhat different on the coarser lattices, particularly at /3 = 6.0, 
but become consistent at P = 6.2. We recall from Fig. 2 of Ref. |j22] that the xLF data at 
P = 6.0 happen to be statistically above their fitted extrapolation, so we see no essential 
disagreement among any of the data sets displayed in our Fig. Q 

The difference between the squared masses of charged and neutral kaons is plotted directly 
in Fig. 121 and is found to be only mildly dependent on (twisted) quark mass over the range 
we are studying. This implies that corrections to Eq. ()39|). arising from higher orders in the 
tmxPT expansion, are small but noticeable. 



15 



Fig. El shows the lattice spacing dependence of the squared mass differences for four mass 
values that span the range of our available data. At leading-order in tm^PT, Eq. fl39p 
indicates that this quantity should be independent of mass, linear in and vanishing in 
the continuum limit. Modulo the unknown higher order effects, Fig. IHl is in reasonable 
agreement with these expectations. In particular, the approximate mass independence is 
evident and the dependence on a? is approximately linear, though a linear fit misses the 
massless prediction at a = by a few (statistical) standard deviations. 

It should be noted from Fig. Elthat even at our smallest lattice spacing, the mass splitting 
of 771^0 — Tnx± ~ 50 MeV is significant relative to the kaon mass itself. However, in terms of 
the difference of mass squared, our results are consistent with the pseudoscalar meson mass 

n n 

splittings in Ref. [22| and compatible with the suggestion of ShindlerjJ] that flavor breaking 
effects in tmLQCD are of a magnitude comparable to "taste" symmetry violations in pseu- 
doscalar meson masses observed with improved staggered fermionsjs^. The appearance of 
sizable lattice spacing effects like this have led some authors to use a power counting scheme 
in which 0{a?) effects are taken to be LO rather than NLoQ,0|, but we will follow Eq. (jZj) 
throughout the present work. 

The decay constant of a charged pseudoscalar meson can be obtained easily from the 



so-called indirect method 



36 



/PS — 2 



(46) 



where the normalization is such that /^r ~ 130 MeV, i.e. larger than the normalization from 
our tmxPT conventions by Unfortunately, the indirect method does not provide easy 
access to the neutral pseudoscalar decay constant, due to mixing with the scalar operator. 
The neutral decay constant is not directly accessible in the laboratory due to the absence 
of flavor-changing neutral currents in the standard model, but it is a quantity that appears 
in the parametrization of some neutral kaon matrix elements (see Ref. for a very recent 
example in the context of tmLQCD). From the point of view of our work, the comparison of 
charged and neutral cases would be able to provide information about how flavor symmetry 
breaking in tmLQCD affects the structure of mesons. Note that if we had chosen a different 
convention in Eq. (^, i.e. interchanging the role of c and s quarks as was done in Ref. [l^ 
which focuses on neutral kaons, then the situation would be reversed: the indirect method 
would have apphed to the neutral kaon and not to the charged kaon. 
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The charged kaon decay constant is plotted in Fig. |3] where we continue to define this 
meson to be the ground state pseudoscalar meson containing one s (anti) quark from the 
heavy doublet and one u (anti) quark from the light doublet. Central values show a hint 
of curvature, but within statistical uncertainties the decay constant is linear in the squared 
meson mass. This is consistent with the tm^PT expression, i.e. the quenched version of 
Eq. ()4U|) where the slope has no logarithmic corrections. (We have verified that a tm^PT 
calculation of the right-hand side of Eq. pUj) also yields Eq. (HO)).) The computations of 
Ref. 3|, using method (ii), also appear in Fig. ^ and are in agreement with the method (i) 
results. 

As is evident from Fig. El there is no visible dependence of the decay constant on lattice 
spacing. For the kaon, fitting the data from three f3 values linearly yields 

roTTiK = 1.25 ^ /i^ = 161 ± 5 MeV . (47) 

Relying on a linear chiral extrapolation, we similarly obtain 

= 142 ± 4 MeV . (48) 

The ratio, 

^ = 1.136(7) , (49) 

Jtt 



agrees nicely with the quenched results from Ref. 
are somewhat larger. Using the lattice spacings derived from Ref. 
us into closer agreement. 



3l| , though the individual decay constants 
3^ or [3^ would bring 



There is also a direct method 



37[ for obtaining the decay constant, though it requires 
input of a renormalization factor for the twisted vector current. Here, we will use the ratio 
of results from the direct and indirect methods to determine this renormalization factor. 
Figure ini shows that the renormalization factor is essentially mass-independent and that it 
becomes closer to unity as a — > 0. The numerical values are comparable to those obtained 
by the authors of Ref. ||2S], and those authors also note that Zy is further from unity in 
tmLQCD than in both standard and boosted lattice perturbation theory. (See their Table 
7.) 

Vector meson masses, referred to here as K* masses since the strange (anti) quark from 
the heavy doublet is combined with a m or (anti) quark from the light doublet, are shown 
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in Fig. [71 as computed from local operators of the form Yll.^i '(p'Jk'^- Within the statisti- 
cal uncertainties, no mass splitting is visible between charged and neutral vector mesons. 
Comparison with data from the xLF collaboration reveals that methods (i) and (ii) lead to 
different K* masses on coarser lattices, and that the distinction vanishes as a 0. The 
sizable uncertainties make chiral extrapolations difficult, particularly for charged mesons. 
Linear fits to the neutral meson masses at each f3 are displayed in Fig. [7| 

It is noteworthy that the neutral K* masses are more precise than the charged K*, just 
as the charged pseudoscalar masses are more precise than the neutral pseudoscalar. In both 
cases, the better precision comes in the channel where the interpolating fields are invariant 
under twisting. Moreover, the absence of large cutoff effects for charged pseudoscalars has 
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been associated with the existence of an exact lattice axial Ward-Takahashi identity. 

Scaling of the neutral vector meson mass with is shown in Fig. |H| The nonvanishing 
dependence on is barely significant with respect to the uncertainties. A linear fit to 
the romps = 1.25 data produces 

rriK' = 970 ± 20 MeV , (50) 

and a linear fit to the linear chiral extrapolations from Fig. [3 yields 

mp = 916 ±20 MeV. (51) 

These quenched values lie above the physical values, but using the lattice spacings derived 
from Ref. or would bring us closer to experiment. 



VI. SCALAR MESON MASSES AND MIXINGS 

Our chosen definition of maximal twist, Eq. ()32|1 . tunes the mixing of vector and axial 
currents for charged mesons, or more precisely, for mesons built from a quark and antiquark 
having twist angles of opposite sign. Charged scalar and pseudoscalar densities do not mix. 
Conversely, there is mixing of neutral scalar and pseudoscalar densities, while neutral vector 
and axial currents do not mix. Our charged vector-axial tuning to maximal twist can differ 
from a neutral scalar-pseudoscalar definition, for example by differing discretization effects. 

Figure IHta) shows four correlators at /3 = 6.2 with our heaviest quark mass: charged 
and neutral scalar and pseudoscalar two-point correlators. The charged correlators cannot 
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mix, and we see a clear ground state exponential behavior for the pseudoscalar and scalar 
mesons, with no contamination between them. The neutral pseudoscalar correlator also 
provides a clear ground state, where the slope in this log plot is slightly steeper than the 
charged expected since we have already established that the neutral meson is heavier 

than the charged meson. The neutral scalar correlator is noticeably different: it maintains 
surprisingly small error bars even far from the source, it displays a kink (change of slope 
on the log plot) near timesteps 12 and 46, and it mirrors the neutral pseudoscalar curve 
between these timesteps. Apparently the quickly-decaying scalar signal is being overcome 
by the pseudoscalar further from the source, i.e. we are seeing scalar-pseudoscalar mixing. 

Figure infb) shows the same four correlators but now with our lightest quark mass. The 
effects are now more dramatic, and a new phenomenon is also observed. The charged scalar 
has a brief signal for the scalar meson near the source, then the correlator becomes negative. 
The neutral scalar similarly has a brief signal, then makes a curious waving shape on the 
graph. To understand this, see Fig. ^1 where the data from Fig. Efb) are replotted on a 
linear scale. The negative contribution to the charged scalar correlator is from the two 
petiole state - ,™,.hed and a. d.scussed Re.. HQ. TKe neut.a. .eala. 

has this two particle state as well, but the correlator is deformed because it apparently also 
has a mixing with the neutral pseudoscalar. 

One could imagine removing the scalar-pseudoscalar mixing by tuning to maximal twist 
directly in this sector, called method (iii) in the notation of Ref. |ll|, but mixing could 
then arise between vector and axial currents. An interesting observation, sketched in Fig. 1 
of Ref. |ll|, is that method (iii) is distinct from method (i) but identical to method (ii) 
up to 0{a?) corrections. Could method (ii) be optimal for both the vector-axial and the 
scalar-pseudoscalar sectors? 

We have computed 100 quark propagators using method (ii) aX [3 = 6.2 and /i^ o = 
0.003608, allowing a direct comparison to our results from method (i). The method (ii) 
value of normal quark mass, m^^o = —0.741546, was obtained from Table 3 in Ref. (27l |. 
Our findings are displayed in Fig. and we see that the scalar-pseudoscalar mixing is still 
apparent, confirming the presence of 0{a?) effects. 
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VII. UNTWISTED STRANGE WITH TWISTED UP AND DOWN QUARKS 

If the twist angle, uj^, of the heavy doublet is set to zero, then the strange quark becomes 
a standard Wilson fermion and its partner can be erased from the action. Exceptional con- 
figurations are typically not a problem for strange quarks, and 0(a) improvement could be 
accomplished via a Sheikholeslami-Wohlert term if desired, though we will use the unim- 
proved Wilson action here. The (m, d) doublet will be kept at maximal twist. 

With this action, the and mesons are exactly degenerate configuration by config- 
uration since one correlator is the Hermitian conjugate of the other. The same is true for 
K~ and . Furthermore, these two pairs are numerically degenerate in the configuration 
average. This can easily be seen in the free quark limit, since the twisted mass in one prop- 
agator cannot contribute to the correlator if the other propagator is Wilson, due to the odd 
number of 75's. Similarly, our tm^PT expression for the mass difference, Eq. fj39|) . explicitly 
vanishes when ujh = ^■ 

Numerical results at /3 = 6.0 for kaon masses obtained with a Wilson strange quark are 
compared to results with a maximally-twisted strange quark in Fig. El For this plot, the 
twisted strange quark is held fixed at aiih,Q = 0.030 and the Wilson strange quark's hopping 
parameter, k, is tuned such that the pseudoscalar mass (obtained from two Wilson propaga- 
tors without any twisting) becomes numerically equal to the charged twisted pseudoscalar 
mass, amps = 0.332(1). The resulting hopping parameter is k = 0.1545 or equivalently 
^?t,o = —0.7634. In both cases, the light quark takes on all four (mj^O;/^«,o) values from 
Table HI 

Fig. El shows that the Wilson strange quark leads to kaon masses that are numerically 
between the charged and neutral kaons with a twisted strange quark. All curves are visibly 
linear in afii^o, though the Wilson strange quark theory has a smaller slope than the twisted 
strange quark theory. Recall that method (i), used here, itself has a smaller slope than 
method (ii). 

Untwisting the strange quark has other effects besides eliminating the mass splitting. 
When an untwisted quark field is combined with a maximally twisted one, parity viola- 
tion induces parity mixing in all (charged and neutral) channels (this can be inferred from 
Eq (jUj)). This is unlike the situation with complete maximal twisting where in some chan- 
nels the parity violation in the action interchanges parity of the correlator but does not mix 
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it. With an untwisted strange quark the extraction of the decay constant becomes a much 
more difficuh problem which is beyond the scope of the present work. 

VIII. SUMMARY 

Twisted mass lattice QCD is a practical method for numerical simulations involving 
light quarks. It has no exceptional configurations, automatic 0{a) improvement, and a 
corresponding version of chiral perturbation theory. However, there are issues of parity and 
flavor symmetry violation effects at non-zero lattice spacing that have to be understood and 
dealt with. 

Quarks come in pairs in tmLQCD, so the best way to implement three-flavor simulations 
requires thought and exploration. In this work, we have considered two doublets at maximal 
twist, where in our quenched simulations the fourth quark is benign. This is in line with the 
two-doublet tmLQCD proposed in Ref. Q|. The chiral perturbation theory is formulated as 
a natural generalization of the existing two-flavor formulation, and used to obtain analytic 
expressions for masses and decay constants. Numerical tmLQCD results for m^^, f^, fn^, fx, 
rrip and mx* are obtained from four twisted quark masses at each of three lattice spacings, 
and are comparable to previous quenched studies with other actions. 

Though dynamical simulations were not performed in this study, that is certainly an 
ultimate goal for QCD phenomenology. Dynamical simulations of the theory with two 
twisted doublets would mean the fourth quark is no longer benign. Identifying it with the 
physical charm quark requires the introduction of a mass splitting within the heavy doublet. 
Progress toward two-doublet dynamical twisted mass simulations is reviewed in Ref. 18|. 
Alternatively, one could avoid an active charm quark by using a mixed action formalism 20| , 
for example with twisted {u,d) and untwisted s quarks in the sea (recall Sec. IVII|1 . and 
Eq. (j2I) used for the valence quarks (so the strange quark's partner is again benign). 

One of the significant twist artifacts found in this work is the mass difference between 
charged and neutral kaons, which vanishes in the continuum limit but remains sizable at 
the lattice spacings studied here, 0.068 fm < a < 0.123 fm. This splitting depends upon 
the particular action that has been chosen; it may be different in other variants of tmLQCD 
or in nonquenched simulations. For example, if only the up and down quarks, not the 
strange quark, are twisted, then this large splitting vanishes however at the price of a 
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more complicated pattern of parity mixing in the correlators. 0{a) errors also arise in that 
scenario, though these could be removed by the addition of a suitable clover operator. 

Another artifact of twisting is the mixing of scalar and pseudoscalar operators when 
the standard definitions of maximal twist are employed. For sufficiently light quarks in the 
quenched approximation, this can be studied through the appearance of negative correlators 
that correspond to the opening of a quenched t]'K channel. 

Notwithstanding the existence of twisted lattice artifacts, we see value in the general 
approach of tmLQCD for applications involving u, d and s quarks. There are a number of 
options for constructing the action including strange quarks, and with systematic studies 
such as this one exploring them we can be hopeful that an optimal approach will be found. 
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APPENDIX A: CURRENTS AND DENSITIES IN THE TWISTED BASIS 

With the generators of SU{2) replaced by those of S'f/(4), the currents and densities in 
the twisted basis of the two-doublet theory are defined in the same way as for the one-doublet 

Q 

theory[7|. At LO, the currents and densities have the same form, mutatis mutandis, as in 
the SU{2) theory. In the physical basis, i.e. in terms of the physical variable Sp/i, they take 
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the forms (with a "hat" denoting the physical basis quantity), 



+ '2(CLm - Ch 



V3 



LO 



pl5 



Slo = cos{J:^u;m)Slo - ^ sin(S"^^)P£o , Plo = cos(S'^o;^)P^o - z sin(S"o;^)52o , 
a, 6Gi^\D, {1,...,15}, D = {3, 8, 15}, (Al) 

where Cp^m = cosUp^m, Sp^m = sintUp^m, P = h, and we use the notation defined in Eqs. (j2Hl- 
OUj) . Note that in the Sf/ (4) theory, P°q and S^^, G ^, do not vanish identically in contrast 
to the SU{2) theory y. 

At NLO, the vector and axial currents are given by 

= V;':^o(l + C) + i^i.2,3,9terms 



+ L. -Tr 



+ W5 -Tr 



[A,, S]_ 9^St - d,J: [Afe, S]_ ) (x'S^ + s/ 



[Afc, S]_ a^St - d,^ [Afc, (ist + sit) 



8aWn 



K = <Lo (1 + C) + ^W,o d.Pto + i.1,2,3,9 terms 



U ^Tr 



- -Tr 



[Afe, S]+ 9^St - [Afc, S]+ j (x'St + Sx't) 
[Afc, S]+ 9^St - [A,, (ist + sit) 



(A2) 



where /c e ^, and 



^ Tr[x'tS + Stxl + ^Tr(itS + Sti) 



(A3) 



We do not give the form of the -^1^2,3,9 terms since each has the same form as in the 
continuum SU{2) theory 41 1. 

Dropping terms proportional to the scalar and pseudoscalar sources, which give rise only 
to contact terms in correlation functions, the scalar and pseudoscalar densities at NLO are 
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given by 



-L8 2SoIV[(AfeEt + EAfe)(x'Et + Ex'^)] - W^g So'IV[(AfeEt + EAfe)(iEt + Ei^)] , 
P'' = P^o (1 + Pi) + + L, BoTr [d^^d^^^A^^^ - EA^)] 

- Ls 2SoTr [(A,Et - EA,)(x'Et + Ex'^)] - Ws i^oTr [(A,Et - EA,)(iEt + eA^)] , 
+ AiH2B^^TriAkH), (A4) 



where A; G ^, and 



V2 



Trp,ED,Et] + ^ Tr(x'tE + Etx) + ^Tr(its + Eti) , 



^Tr(x'tE-Etx')-^Tr(itS-Eti) 



(A5) 



To write the NLO currents and densities in the physical basis, we need the results 
Tr(DtED^Et) ^ Tr(Dt5^,,D^E;j , 

IV(x'^E±Etx) = 



4M' go AAM' ( 1 q8 i 1 Cl5 \ /' i ^Jo-ti'I 



4M' 



4AM'/ 1 f)8 



p8 J L Pl5 



+ 0(M;ep), 



signj 



Tr(i^E ± E^i) 



^Slo (+sign) 



+ 0{a€p) , p^l,h, 



where 



(A6) 



(A7) 



M' = (m; + MO/2 , AM' = m; - m;; , 

which allow us to express C, Vi, and Vf in terms of the physical fields. 

Next we write the L5, Lg, W5, Ws, and terms in the physical basis. For the L5 terms, 
we need the results 



± -Tr 
2 



([A„, E]^ a^Et - a^E [A„, E]^) (x'Et + Ex'^)" 

cos(A"cj„) - r]ab sin(A"w„)l^g (upper sign) 
cos(A"(x;^)i25 - Vab sin( A"u;^) (lower sign) 

cos(S"cj„)^t -«sin(S"cj„,)P£, 



IV[9^E9^Et(A„Et±EA„)] 



cos(E»u;„)P£^ - isin(E»u;^)^£^ 



sign) 
sign) 



(A8) 
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Atlii SmA Sm 



where a, b & -^V^), and 

= ^ (Aa ( xLS^ + ^ + (S ^ St, 

^2, = - (Aa(5^Sp,St,a^Sp, + h.c.)) , Pl^ = (A,(9^Sp,St,9^S,, - h.c.)) , (A9) 
and for the La terms we need the results 



Tr[(A.St±SA,)(x'St + Sx't 



Sls = - {Ki^phX'L^ph + h.c.) 



cos(S"cu„)5£^ - I sin(S'^c^^)P£^ (+ sign) 
cos(S"a;„)P£^ - z sin(S»^^)5£^ (- sign) 

Pl^ = (A,(Sp,x;lSp, - h.c.)) , a G . 

(AlO) 



By replacing x' with A, and L5 g with Ws^s, the and Ws terms can be expressed in the 
physical basis using the same results above for the L5 and Lg terms. 

Lastly, the H2 term contributes only in the flavor-diagonal case, i.e. when the flavor 
index A; = 3, 8, 15. Since we will not be using the flavor- diagonal currents and densities, we 
do not give results for the flavor- diagonal cases here. 

To conclude this appendix, we provide the explicit expression for the axial current in the 
physics basis at NLO, using the twist angles determined in Sec. 1111 Bl 

8Woa, 



Aa _ Aa 



-WiQ cos(SVm) ?,m.{/VuJm)r]ahV ^qO- + ^) + ^1,2,3,9 terms 



+ ^^W^io cos(A'^o^^) cos(S'^u;^)9^P^ - i sin(S"a;^)9^S2o 



A 



U - ^^W^io cos(S'^a;™) sin(A'^cu^)r/,,y2^ 



+ [L5 ^ W,] . (All) 



The term C is the same as in Eq. ()A3|1 . but now given in the physical basis. 
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TABLE I: The parameters used for simulations in this work. Lattice spacings are taken from 
Ref. J^]. Each {amp^Q^a^p^o) pair is the result of tuning to maximal twist with method (i) as 
discussed in Sec. IIVI The subscript p = l,h is used in the text to distinguish the "light" quark 
doublet from the "heavy" quark doublet, but for purposes of numerical tuning in this table there 
is no distinction. The twist angle was obtained from Eq. (|41jl . 



p 


a [fm] 


#sites 


^configurations 


anip^o 




twist angle (degrees) 


5.85 


0.123 


20^ X 40 


600 


-0.8965 


0.0376 


90.0±0.3 










-0.9071 


0.0188 


90.2±0.6 










-0.9110 


0.01252 


90.6±0.8 










-0.9150 


0.00627 


90.6±1.6 


6.0 


0.093 


20^ X 48 


600 


-0.8110 


0.030 


90.4±0.4 










-0.8170 


0.015 


91.0±0.7 










-0.8195 


0.010 


92.5±1.0 










-0.8210 


0.005 


95.5±2.1 


6.2 


0.068 


28^ X 56 


200 


-0.7337 


0.021649 


89.1±0.8 










-0.7367 


0.010825 


87.3±1.8 










-0.7378 


0.007216 


86.3±2.8 










-0.7389 


0.003608 


86.4±4.5 
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TABLE II: Numerical values from our simulations. These are also shown graphically in the fig- 
ures. The rows with superscripts a and b refer to (a/U/j^0)O^W,o)=(0-015, 0.005) and (0.010,0.010) 
respectively. 

/? afii^o + aHh,o (amps)^ a/ps Zy amy 



charged neutral 



charged neutral 



5.85 0.0752 
0.0564 
0.05012 
0.04387 
0.0376 
0.03132 
0.02507 
0.02504 
0.01879 
0.01254 



0.1841(4 
0.1388(4 
0.1236(4 
0.1085(4 
0.0937(3 
0.0784(10) 
0.0633(3 
0.0633(3 
0.0484(2 
0.0327(2 



0.2262(12) 
0.1827(14) 
0.1683(15) 
0.1536(19) 
0.1402(15) 
0.1259(16) 
0.1112(19) 
0.1117(17) 
0.0969(22) 
0.0824(33) 



0.1127(8) 

0.1064(8) 

0.1041(9) 

0.1019(9) 

0.0996(9) 

0.0972(13) 

0.0947(9) 

0.0946(9) 

0.0921(9) 

0.0892(10) 



0.620(4) 

0.615(4) 

0.614(4) 

0.613(4) 

0.607(4) 

0.602(10) 

0.602(5) 

0.601(5) 

0.601(5) 

0.597(7) 



0.625(5) 

0.591(8) 

0.580(10) 

0.573(13) 

0.553(14) 

0.539(18) 

0.523(26) 

0.522(25) 

0.497(37) 

0.456(58) 



0.622(3) 

0.591(5) 

0.582(6) 

0.575(8) 

0.563(7) 

0.555(9) 

0.548(10) 

0.547(10) 

0.539(12) 

0.527(18) 



6.0 0.060 


0.1106(4) 


0.1260(6) 


0.0858(7) 


0.661(4) 


0.488(4) 


0.484(3) 


0.045 


0.0829(4) 


0.0986(7) 


0.0811(7) 


0.656(5) 


0.463(6) 


0.462(4) 


0.040 


0.0738(4) 


0.0898(8) 


0.0796(7) 


0.654(5) 


0.455(7) 


0.456(5) 


0.035 


0.0646(5) 


0.0801(9) 


0.0780(8) 


0.654(6) 


0.444(10) 


0.450(7) 


0.030 


0.0558(4) 


0.0724(7) 


0.0762(7) 


0.649(5) 


0.437(8) 


0.441(6) 


0.025 


0.0468(4) 


0.0641(8) 


0.0745(8) 


0.647(6) 


0.429(10) 


0.437(7) 


0.020" 


0.0378(4) 


0.0550(10) 


0.0726(8) 


0.646(7) 


0.418(15) 


0.432(9) 


0.020'' 


0.0378(4) 


0.0559(9) 


0.0727(8) 


0.644(6) 


0.421(14) 


0.433(9) 


0.015 


0.0290(3) 


0.0471(12) 


0.0706(8) 


0.644(7) 


0.412(19) 


0.430(11) 


0.010 


0.0198(3) 


0.0383(18) 


0.0680(9) 


0.637(11) 0.407(27) 


0.426(15) 


6.2 0.043298 


0.0585(4) 


0.0640(6) 


0.0614(7) 


0.692(10) 


0.362(4) 


0.360(3) 


0.032474 


0.0441(4) 


0.0497(6) 


0.0582(8) 


0.689(11) 


0.345(6) 


0.344(5) 


0.028865 


0.0393(4) 


0.0451(7) 


0.0571(8) 


0.689(11) 


0.340(8) 


0.340(5) 


0.025257 


0.0346(5) 


0.0406(7) 


0.0562(8) 


0.688(11) 


0.335(9) 


0.337(6) 


0.02165 


0.0298(4) 


0.0358(7) 


0.0547(8) 


0.686(12) 


0.328(9) 


0.329(6) 


0.018041 


0.0250(4) 


0.0313(7) 


0.0536(9) 


0.685(13) 


0.322(11) 


0.325(7) 


0.014433 


0.0203(4) 


0.0269(8) 


0.0525(10) 0.684(13) 


0.317(14) 0.324(8) 


0.014432 


0.0203(4) 


0.0268(8) 


0.0523(10) 0.684(14) 


0.317(14) 0.322(8) 
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FIG. 1: Pseudoscalar meson mass squared as a function of the sum of quark and antiquark twisted 



mass parameters. Subscripts I and h indicate the hght and heavy doublets. Result^ 
method (i) are from the present work; results labelled by method (ii) are from Ref. 



abelled by 



27l | and have 



equal masses for the quark and anti-quark. Straight lines are linear fits to the data from method 
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FIG. 2: The difference between charged and neutral squared pseudoscalar meson masses as a 
function of the sum of quark and antiquark twisted masses. Subscripts I and h indicate the Ught 
and heavy doublets. 
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FIG. 3: The difference between charged and neutral squared pseudoscalar meson masses as a 
function of squared lattice spacing, for selected values of the charged meson mass. 
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FIG. 4: The pseudoscalar meson decay constant as a function of the squared charged pseudoscalar 
meson mass. Results labelled by method (i) are from the present work; results labelled by method 
(ii) are from Ref. [28]. Straight lines are linear fits to the data from method (i). 
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FIG. 5: Scaling of the pseudoscalar decay constant for four choices of the quark mass. 
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FIG. 6: The renormahzation factor associated with the pseudoscalar meson decay constant. 
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FIG. 7: Vector meson mass as a function of the squared charged pseudoscalar meson mass. Results 
labelled by method (i) are from the present work; results labelled by method (ii) are from Ref . 25 1 . 
Straight lines are linear fits to the neutral data from method (i). 
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FIG. 8: Scaling of the neutral vector meson mass for four choices of the quark mass. The straight 
line is a linear fit to the data having romps = 1.25. 
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FIG. 9: Scalar and pseudoscalar correlation functions for our (a) heaviest and (b) lightest quarks 
at /? = 6.2. The notation 75 and 1 refers to the physical basis, defined using method (i). 
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FIG. 10: The data from Fig. ini(b), replotted on a linear scale. 



0.0035 p 
0.003 r 
0.0025 r 
0.002 r 
0.0015 r 
0.001 r 
0.0005 r 
- 
-0.0005 r 

-o.ooiL 



IE 
IE 



jy charged, method (i) 
y^, charged, method (ii) 
jy neutral, method (i) 
y^, neutral, method (ii) 
1, neutral, method (i) 
1, neutral, method (ii) 
1, charged, method (i) 
1, charged, method (ii) 




60 



FIG. 11: The data from Fig. 1101 compared to computations using method (ii). 
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FIG. 12: Squared pseudoscalar meson masses with one strange quark/antiquark and one light 
quark/antiquark, plotted as a function of the light quark's twisted mass. 



